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ABSTRACT 

We predict the fraction of dark halo lenses, that is, the fraction of lens systems produced by the gravitational 
potential of dark halos, on the basis of a simple parametric model of baryonic compression. The fraction of dark 
halo lenses primarily contains information on the effect of baryonic compression and the density profile of dark 
halos, and is expected to be insensitive to cosmological parameters and source population. The model we adopt 
comprises the galaxy formation probability po{M) which describes the global efficiency of baryonic compression 
and the ratio of circular velocities of galaxies to virial velocities of dark halos 7,, = Vc/vvk which means how the 
inner structure of dark halos is modified due to baryonic compression. The model parameters are constrained 
from the velocity function of galaxies and the distribution of image separations in gravitational lensing, although 
the degeneracy between model parameters still remains. We show that the fraction of dark halo lenses depends 
strongly on 7,, and the density profile of dark halos such as inner slope a. This means that the observation of 
the fraction of dark halos can break the degeneracy between model parameters if the density profile of dark halo 
lenses is fully settled. On the other hand, by restricting 7,, to physically plausible range we can predict the lower 
limit of the fraction of dark halo lenses on the basis of our model. Our result indicates that steeper inner cusps of 
dark halos ia>l .5) or too centrally concentrated dark halos are inconsistent with the lack of dark halo lenses in 
observations. 

Subject headings: cosmology: theory — dark matter — galaxies: formation — galaxies: halos — galaxies: 
clusters: general — gravitational lensing 



1. INTRODUCTION 

Strong gravitational lensing offers a powerful probe of 
the matter distribution in the universe. So far ~ 60 lensed 
quasars are known, and their properties are summarized by 
CfA/ Arizona Space Telescope Lens Survey (CASTLES'). 
Most of these lens systems have image separations 9 < 3". The 
cold dark matter (CDM) scenario, however, predicts sufficiently 
cuspy dark halos (e.g., Navarro, Frenk, & White 1996, 1997), 
thus dark halos can produce the significant amount of multi- 
ple images even at 6* > 3" (Wyithe, Turner, & Spergel 2001; 
Keeton & Madau 2001; Takahashi & Chiba 2001; Li & Os- 
triker 2002; Oguri et al. 2002). It has been unclear whether 
the distribution of image separations should be computed based 
on galaxies (using the luminosity function and the density pro- 
file of galaxies; Turner, Ostiiker, & Gott 1984; Turner 1990; 
Fukugita & Turner 1991; Fukugita et al. 1992; Maoz & Rix 
1993; Kochanek 1996; Chiba & Yoshii 1999) or dark halos (us- 
ing the mass function and the density profile of dark halos; 
Narayan & White 1988; Cen et al. 1994; Wambsganss et al. 
1995; Kochanek 1995; Maoz et al. 1997; Wambsganss, Cen, 
& Ostriker 1998; Mortlock & Webster 2000; Wyithe, Turner, 
& Spergel 2001; Keeton & Madau 2001; Takahashi & Chiba 
2001; Oguri et al. 2002). The distribution of image separations, 
however, clearly indicates that the use of only dark matter prop- 
erties cannot match observations. That is, the observed distribu- 
tion of image separations is never reproduced from the theoret- 
ical calculation using the mass function of dark halos and only 
one population for lensing objects (e.g., Li & Ostriker 2002). 
The possible solution to explain the distribution of image sep- 
arations is the modification of inner structure of dark halos by 
introducing baryonic cooling (Keeton 1998; Porciani & Madau 



2000; Kochanek & White 2001; Keeton 2001; Sai'bu, Rusin, & 
Ma 2001; Li & Ostriker 2002). In this picture, inside low mass 
halos baryons are efficiently compressed and form sufficiently 
concentrated galaxies, while inside larger halos such as group- 
or cluster-mass halos global baryon cooling hardly occurs and 
thus the inner structure of dark halos remains unmodified. 

Then a question comes to our mind: what is the fraction of 
dark halo lenses? Here the term "dark halo lenses" is used to 
describe lenses which are produced by the gravitational poten- 
tial of dark halos. In other words, the lens objects of dark halo 
lenses are dark halos which have no central galaxies or have 
central galaxies but they are too small to dominate in gravi- 
tational lensing. Dark halo lenses exhibit characteristic prop- 
erties such as the small flux ratios and the detectable odd im- 
ages (Rusin 2002), thus can be distinguished from usual galaxy 
lenses. Since dark halos are expected to have a steep cen- 
tral cusp, the significant amount of dark halo lenses should 
be observed. But so far no confirmed dark halo lens system 
is observed in strong gravitational lensing survey. The excep- 
tion is arc statistics in rich clusters (Bartelmann et al. 1998; 
WilHams, Navarro, & Bartelmann 1999; Meneghetti et al. 2001 ; 
Molikawa & Hattori 2001; Oguri, Taruya, & Suto 2001; Oguri 
2002), but the known cluster lenses were all found by search- 
ing for lenses in detail after identifying a rich cluster In the 
surveys which first identify source objects and see whether they 
are lensed or not, it seems that dark halo lenses have not been 
observed yet: statistical argument (Kochanek, Falco, & Munoz 
1999) and individual properties (Rusin 2002) imply that current 
ambiguous quasar pairs are likely to be binary quasars. Even 
known lensing systems in clusters, such as Q0957H-561 (Walsh, 
Carswell, & Weymann 1979), are produced mainly by a galaxy 
in the cluster. The cluster potential contributes to lensing only 
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as a perturbation. Therefore, it should be checked whether the 
lack of dark halo lenses in observations really reconciles with 
the theoretical prediction. 

Description of baryonic effects needs detailed models for 
the star formation and feedback (e.g., Cole et al. 2000). In- 
stead, in this paper we predict the expected fraction of dark 
halo lenses on the basis of a simple (minimal) parametric model 
(Kochanek 2001). This model comprises the formation proba- 
bility of galaxies, Pg(M), and the ratio of circular velocities of 
galaxies to virial velocities of dark halos, 7,.. The former de- 
scribes the global efficiency of baryonic compression, and the 
latter models the modification of inner structure of dark halos 
due to baryonic compression. The model parameters are cho- 
sen so as to reproduce the velocity function of galaxies and the 
distribution of image separations. Although there remains the 
strong degeneracy between model parameters, this degeneracy 
can be broken from the observation of the fraction of dark halos 
if we fix the density profile of dark halos. On the other hand, by 
restricting a range of 7y from various theories and observations, 
we can also derive the lower limit of the fraction of dark halo 
lenses. Our main finding is that steep iimer slopes of dark halos 
(a > 1 .5) or too centrally concentrated dark halos are incon- 
sistent with the lack of dark halo lenses in observations, even 
if various uncertainties are taken into account. Although this 
constraint on the density profile is not so severe, a large lens 
sample obtained by e.g., Sloan Digital Sky Survey (SDSS) can 
put tighter constraints on the density profile of dark halos as 
well as the model of baryonic compression. 

The plan of this paper is as follows. In §2, we de- 
scribe the model of baryonic compression. Section 3 is de- 
voted to constrain the model parameters, and §4 presents 
our predictions for the fraction of dark halo lenses. Finally, 
we surmnarize conclusions in §5. Throughout this paper, 
we assume the lambda-dominated cosmology i^lo,^o,h,a^) = 
(0.3,0.7,0.7,1.04), where the Hubble constant in units of 
100kms"'Mpc"' is denoted by h. As shown below, however, 
our results are quite insensitive to a particular choice of cosmo- 
logical parameters. 

2. THEORETICAL MODEL 

2.1. Effects ofBaryons on Dark Halos 

The obvious difference between velocity functions of galax- 
ies and dark halos at high velocity (e.g., Gonzalez et al. 2000) 
indicates that the efficiency of the baryon compression changes 
from galaxy-mass scale to group- and cluster-mass scales. If the 
mass of dark halos is sufficiently large, the baryon cooling time 
Tcooi which increases with halo mass (e.g.. Cole et al. 2000) be- 
comes larger than the age of dark halos which slightly deceases 
with halo mass (e.g., Lacey & Cole 1993; Kitayama & Suto 
1996). Thus baryons inside such a massive dark halo do not 
form so large galaxies as to collect most of baryons in the dark 
halo. This yields a steep cutoff at high velocity in the galaxy ve- 
locity function. To model this, following Kochanek (2001), we 
introduce the probabihty Pg(M) that a sufficiently large galaxy 
which collects most of internal baryons is formed inside a halo 
of mass M. For example, galaxy-mass halos with mass M 
usually have large galaxies with mass Mgai ~ i^b/Clo)M. On 
the other hand, groups or clusters also have galaxies but their 
mass are small even for galaxies which lie at the center of ha- 
los, Mgai < iQ.b/^o)M, here in this case M > lO'^/i^'M© (e.g., 
Yoshikawa et al. 2001). Therefore we regard galaxies in groups 
or clusters as "substructures" and do not take into account. We 



adopt a parametric model (Kochanek 2001) 



Pg(M) = <^ 
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(M < Mh), 
(M > Mh). 
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Although the overall factor of Pg(M) should become another pa- 
rameter, we neglect it. The effect of the overall factor of Pg{M) 
on the fraction of dark halo lenses will be discussed in §4. 

Also, there is a difference in the slopes of velocity functions 
of dark halos and galaxies at small velocity (e.g., Gonzalez et 
al. 2000; Nagamine et al. 2001), and this demands the modi- 
fication of Pg{M) at low M. We, however, neglect this effect 
because we are not interested in low velocity galaxies which 
hardly contribute to lensing statistics. 

The baryon compression also changes inner structure of dark 
halos (e.g., Blumenthal et al. 1986; Mo, Mao, & White 1998). 
Before the baryon cooling occurs, we assume that the density 
profile of dark halos is well described by the one-parameter 
family of the form (Zhao 1996; Jing & Suto 2000) 



Pir) = 



3-a ' 



(2) 



where = rvir/cvk and Cvir is the concentration parameter. We 
adopt the mass and redshift dependence reported by Bullock et 
al. (2001): 



Cyir(M,z): 



M 



-0.13 



(3) 



for a = I, and generalize it to a ^ 1 by multiplying the fac- 
tor (2 -a) (Keeton & Madau 2001). We also take account of 
scatter of the concentration parameter which has a log-normal 
distribution with the dispersion of (Tc = 0. 1 8 (Jing 2000; Bullock 
et al. 2001). The characteristic density 6^ can be computed us- 
ing the spherical collapse model (see Oguri et al. 2001). While 
the correct value of a is still unclear, the existence of a cups 
with i ^ a ^1-5 has been established in recent N-body simu- 
lations (Navarro et al. 1996, 1997; Moore et al. 1999; Ghigna 
et al. 2000; Jing & Suto 2000; Klypin et al. 2001; Fukushige & 
Makino 2001, 2002; Power et al. 2002). 

Next consider the modification of inner structure of dark 
halos induced by baryonic compression. In general, cooled 
baryons (galaxies) are more centrally concentrated than dark 
matters. Although complicated physical models are needed to 
know the modified mass distribution, we simply assume that the 
mass distribution of galaxies is well described by the Singular 
Isothermal Sphere (SIS) approximation: 



p(r) = 



(4) 



where ct is a one-dimensional velocity dispersion and is related 
to the circular velocity as Vc = \/2(t. We also assume that Vc 
as a function of halo mass M can be described as 



7v = 



Vvir(M)' 



(5) 



where Vvii(M) is the halo virial velocity defined by Vvii(M) = 
^/GMjr\^ and 7,, is an arbitrary constant. The SIS approxima- 
tion for galaxies is consistent with several observations such as 
dynamics (e.g., Rix et al. 1997) and gravitational lensing (e.g., 
Rusin & Ma 2001; Cohn et al. 2001). 



FRACTION OF DARK HALO LENSES 



3 



2.2. Velocity Function of Galaxies 

Once the probability Pg(M) and the galaxy circular veloc- 
ity Vc(M) are given, we can calculate the velocity function of 
galaxies from the mass function of dark halos: 



dn 
dvc 



(Ve(M))=Pg(M) 



dM 



dn 
dM 



(M). 



(6) 



For the mass function dn/dM, we adopt a fitting form derived 
by Sheth & Tormen (1999): 



dn, 



ST 



dM 
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1 + 
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0.707 Sliz) J 



1.414 



, Po So(z) 
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da, 
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dM 



exp 



0.7075l{z} 



(7) 



This fitting form coincides more accurately with numerical 
simulations than the analytic mass function derived by Press 
& Schechter (1974). The velocity function is sometimes ex- 
pressed in terms of log Vc as 

dn 

^(yc) = - . (8) 

c/logVc 

In this calculation we neglect the contribution from "substruc- 
tures" (i.e., galaxies in groups and clusters), because this mainly 
changes the normahzation of the velocity function which we do 
not use as a constraint on the model parameters (see §3.1). Sub- 
structures, however, may affect the fraction of dark halo lenses 
directly, because it changes only the number of galaxy lenses. 
Therefore we consider the effect of substructures in predicting 
the fraction of dark halo lenses (see §4). 

2.3. Lensing Probability Distribution 

Bearing the picture described in §2.1 in mind, we calculate 
the probability of gravitational lensing caused by bright galax- 
ies and dark halos separately: 

Pg^(>^;zs,i-) = 

Pd,A(>0;zs,L) = 

{ 1 -pg(M)} aNFw£:T-(l +Zl)^ :ttt j(10) 

where zs is the source redshift, Zh is the lens redshift, L is the 
luminosity of the source, and dn/dM is the comoving number 
density of dark halos (eq. [7]). Lensing cross sections ctsis and 
C7NFW are given by 

2 

(11) 



JMpg(M)asisB^(l+ZL)'^, (9) 
dzh dM 



Jo JM„ 



(Tnfw 



(12) 



where is the radius of the radial caustic at source plane (e.g., 
Schneider, Ehlers, & Falco 1992) and /Jql, £>os> and Z)ls are 
the angular diameter distances to the lens, to the source, and be- 
tween the lens and source, respectively. The one-dimensional 
velocity dispersion u in equation (11) can be represented as a 
function of halo mass M: 

Vc 7„Vvir(M) 



where equation (5) is used. The lower limit of integral by mass, 
Mmin, is determined by solving the equation = 0(Mjain,zs,ZL)- 
The magnification bias (Turner 1980; Turner et al. 1984) is in- 
cluded in B as follows: 

1 f 1 



B-- 



d'ri^{z5,L/ii{fi))- 



(14) 



0'lens^(ZS,i') . ' multi 

where $(zs,i') is the luminosity function of sources and ^lijj) is 
the magnification factor at -q. 

The total lensing probability with image separation larger 
than 6 is given by 

P{> ^;zs,L) = Pgai(> ^;zs,Z.)+Pdark(> e;zs,L). (15) 

3. CONSTRAINTS ON THE MODEL 

3.1. Velocity Function of Galaxies 

The velocity function is valuable because it is easy to handle 
theoretically, and is useful to test the model of galaxy formation 
and cosmology (Cole & Kaiser 1989; Shimasaku 1993; Gonza- 
lez et al. 2000; Kochanek & White 2001). Using the model de- 
scribed in §2.2, we can calculate the velocity function of galax- 
ies from the mass function of dark halos. By comparing the 
theoretical velocity function with observed velocity functions, 
one can constrain the model parameters such as 5b, Mb, and 7v. 
The definite theoretical prediction, however, needs the correct 
value of cosmological parameters such as ag. Therefore we use 
the normaUzed velocity function 

w(200km/s) 

instead of vE'(vc). Moreover, by using V'(vc) we can neglect the 
effect of the overall factor of Pg{M). For the observed veloc- 
ity functions, we use five velocity functions derived by Gonza- 
lez et al. (2000): velocity functions derived from Southern Sky 
Redshift Survey (SSRS2), Automatic Plate Measuring facil- 
ity survey (APM), United Kingdom Schmidt Telescope survey 
(UKST), Las Campanas Redshift Survey (LCRS), and K-band 
survey by Gardner et al. (1997). We restrict the comparison 
with observations in the region 200km/s < Vc < 500km/s be- 
cause of the following two reasons. One reason is that observed 
velocity functions are not reliable at the high velocity region 
(vc > 500km/s) (Gonzalez et al. 2000). The other reason is the 
obvious difference between observed and theoretical number 
density of galaxies at the low velocity region (vc < 200km/s). 
The reason of this difference is that in our model we neglect 
effects of e.g., supemovae feedback (Dekel & Silk 1986) which 
substantially suppress the number of galaxies in low-mass ha- 
los (e.g., Gonzalez et al. 2000; Nagamine et al. 2001). However 
we are interested in the gravitational lensing with angular sepa- 
ration larger than current angular resolution {0 ~ 0.3"), thus we 
can safely neglect such low- velocity galaxies which have Uttle 
impact on our lensing results. 

Figure 1 shows constraints on model parameters from ob- 
served velocity functions. Contours are calculated from ra- 
tios of combined UkeUhood C oc ne''^P(~2X?)- Errors of ob- 
served velocity functions are estimated from errors of fitting 
parameters of velocity functions (see Gonzalez et al. 2000). 
As easily seen from the upper panel of Figure 1, we can put 
useful constraints on and Vh = Vc(Mh). The best fit param- 
eter set is (^hjV'h) ~ (0.78,238km/s). This means that there 
still remains strong degeneracy between M^ and 7^: parame- 
ter combinations which yield the same Vh cannot be discrimi- 
nated from the velocity functions. This fact is also shown in 
the lower panel of Figure 1. One of our best fit parameter 
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sets (^h,^^h,7v) = (0.78, 1O^^^~'M0, 1.67) is compared with ob- 
served velocity functions in Figure 2. It is clearly shown that 
the prediction of our model shows good coincidence with ob- 
served velocity functions. 
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Fig. 2. — One of our best fit models for the velocity function is compared 
with the five observed velocity functions (Gonzalez et al. 2000). Here normal- 
ized velocity functions (eq. [16]) are compared. Errorbars are estimated from 
the error of fitting parameters of velocity functions. The model parameters are 
chosen so as to satisfy constrains in Figure 1: <5h = 0.78, A/h = lO'^/i-'A/©, 
and 7v = 1.67. 



3.2. Image Separation Distribution 

The distribution of image separations in strong gravitational 
lensing becomes another test of our model. The distribution 
is calculated from equation (15). For the observed lens sam- 



ple, we consider the Cosmic Lens All-Sky Survey (CLASS; 
Helbig 2000). This sample is complete at image separations 
0.3" <e< 15" (Helbig 2000; Phillips et al. 2001), and has 18 
lenses among ~ 12,000 radio sources. Sources have the flux 
distribution dn/dS oc 5"^ ' (Rusin & Tegmark 2001). The mean 
redshift is estimated to be (zs) = 1 -27, although the redshift dis- 
tribution of sources is still poorly understood (Marlow et al. 
2000). As in the case of the velocity function, actually we con- 
sider the normaUzed image separation distribution: 

^^>''^^'^^=P^03WO' ^''^ 
instead of usual probability distribution. This is because the ab- 
solute probability may suffer from uncertainties of source red- 
shifts, magnification bias, and cosmological parameters. On 
the other hand the distribution p(> 6) mainly contains informa- 
tion on the shape of the mass function and effects of baryonic 
compression (Kochanek 2001; Kochanek & White 2001; Kee- 
ton 2001). And the distribution p{> 6) is quite insensitive to 
the source population and cosmological parameters, thus can 
be used for samples in which the source population is unknown 
such as CASTLES sample. 

The distribution of image separations is tested by using a 
Kolmogorov-Smirnov (KS) test. For the observed distribution, 
we use the distribution of CLASS survey. The result is shown 
in Figure 3. We plot contours in Mh-7v plane by fixing 5h = 0.78 
which is constrained from the velocity function (see Figure 1). 
From this figure, we find that the constraint from the distribu- 
tion of image separations is consistent with the constraint from 
the velocity function. The exception is the case of a = 1.5, 
where a is the irmer slope of the density profile of dark halos 
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Fig. 3. — Constraints on model parameters from the distribution of image separations in strong gravitational lensing. Here we have assumed = 0.78. 68% 
(solid), 95% (dotted), and 99% (dashed) contours are calculated using Kolmogorov-Smimov test. 68% constraint from the velocity function (see Fig. 1) is also 
plotted by thick Unes. For the iimer slope of the density profile of dark halos (eq. [2]), we consider both a = 1 .0 (left) and a = 1 .5 (right). 



(eq. [2]). In this case, the lower value of 7v becomes incon- 
sistent with the observation, because the lensing cross section 
for SIS is (Tsis oc 7^ (see eq. [11]) and the lower value of 71, 
means that the contribution of dark halo lenses becomes more 
significant. That is, when 71, is low, the number of dark halo 
lenses is comparable to that of galaxy lenses, thus a steep cutoff 
in the distribution of image separations never appears. There- 
fore, we conclude that the constraint from the distribution of 
image separations is consistent with the constraint from the ve- 
locity function, except for the cutoff of KS probability at low 7v 
which may appear when a is large. We also compare the distri- 
bution of our model with the observed distribution in Figure 4. 
The model parameters are chosen so as to satisfy the constraint 
from the velocity function, and same as those used in Figure 
2. We also plot the distribution assuming that all lens objects 
are well approximated by SIS. This assumption corresponds to 
Pg(M) = 1 for all M. It is obvious that our model well repro- 
duces the observed distribution while the distribution of all SIS 
assumption is far from the observed distribution, as pointed out 
by Keeton (1998). 
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Fig. 4. — The distribution of image separations in gravitational lensing (eq. 
[17]). The model parameters are same as Figure 2. The cases with a = 1.0 
(dotted) and a = 1.5 (dashed) are plotted. The observed distribution in CLASS 
survey is shown by thick solid Une. The thin solid line indicates the theoreti- 
cal distribution assuming all lens objects are SIS lenses. This corresponds to 
Pg(M) = 1 for all M. 



3.3. Various Constraints on the Circular Velocity of Galaxies 

In §3.1 and §3.2, we constrained model parameters from the 
velocity function of galaxies and the distribution of image sep- 
arations. The strong degeneracy between Mh and 7^, however, 
still remains. Therefore in this subsection we try to put con- 
straints on 7v from various theories and observations. We con- 
sider following three constraints. 

First, Seljak (2002a) gave values of 7^ at several halo mass 
derived from the observation and analysis of galaxy-galaxy 
lensing (McKay et al. 2002; Guzik & Seljak 2002) and TuUy- 
Fisher/fundamental plane relations. The observation of galaxy- 
galaxy lensing allows us to determine the mass of dark ha- 
los and its relation to the luminosity of galaxies. TuUy- 
Fisher/fundamental plane relations are used to derive the galaxy 
circular velocity from its luminosity. He found that 7v is signif- 
icantly larger than 1, 7^ ~ 1.8 around L,. He also found the 
decrease of 7v fromL* to 7L*, 7v ~ 1.4 at 7L*. 

Second, Cole et al. (2000) semi-analytically calculated av- 
erage values of 7^ for galaxies with -20 < M/-51og/! < -18. 
They calculated 7,, for various parameter sets of semi-analytic 
model, and found that most of parameter sets predict 7v ~ 1 .3 - 

1 .4. The most important parameter for 7v is the baryon density 
J7b, but even in the extreme cases they examined, fib = 0.01 and 
0.04, 7,, becomes 1.13 and 1.96, respectively. 

Third, Mo et al. (1998) considered the analytic disk forma- 
tion model on the basis of adiabatic compression (e.g., Blu- 
menthal et al. 1986). In their model, 7v depends on the mass 
of cooled baryons and the angular momentum of galactic disc. 
More specifically, their model is characterized by following 
three parameters; the concentration parameter c, the ratio of 
disk mass to halo mass mj, and the angular momentum A' (see 
Mo et al. 1998). We assume that the concentration parameter is 
related to halo mass as in equation (3). Then 7,, at fixed mass de- 
pends on and A'. We examine following two extreme cases: 
(md, A') = (0.02,0.02) and (0.02,0.2). Most of parameter sets 
which are physically reasonable predict the value of 7v between 
these two extreme cases. 



6 



OGURI 



6 I 1 1 — I — I I I 1 1 1 1 1 — I — I I I 1 1 1 1 1 — I — I I I 1 1 

5 - A Seliak(2002a) 

4 - X Cole et al.(2000)- 

Shading Mo et al.(1998) 




0.7 - 



■ 10" 10'2 1013 iQi* 

M[h-iMg] 

Fig. 5. — Constraints on 7v from various models, (i) Seljak (2002a) gave 
values of 7v at several halo mass derived from the observation and analysis of 
galaxy-galaxy lensing (McKay et al. 2002; Guzik & Seljak 2002) and TuUy- 
Fisher/fundamental plane relations. Errorbars indicate Icr level error, (ii) Cole 
et al. (2000) semi-analytically calculated the average values of 7v for various 
parameter sets of semi-analytic model. The errorbar in the direction of 7v-axis 
means that aH parameter sets predict 7v within this range, (iii) In the analytic 
disk formation model of Mo et al. (1998), 7v depends on the mass of cooled 
baryons and the angular momentum of galactic disc. Any choice of these two 
parameters which is physically allowed, however, results in the value of 7v in 
the shading region. 

The result is summarized in Figure 5. Particularly we are in- 
terested in the strong gravitational lensing, thus we focus on the 
range 1O'^/!"'M0 <M< 1O'^/i"'M0 where gravitational lens- 
ing with separations 6 > 0.3" is most efficient. In this range 
of halo mass, various constraints indicate that 7i, is restricted to 
1 < 7v < 2. In particular, the upper bound of jy is quite robust 
because 7v is closely related to the amount of cooled baryons; 
to produce ^ 2, we need the extraordinarily large amount 
of cooled baryon which exceeds the global baryon mass ratio 
(Seljak 2002a). Figure 5 also indicates that 7,, should slightly 
depend on halo mass. The mass dependence of 7v, however, 
is not so important because our interest is restricted in a narrow 
mass range. Therefore we can safely assume that 7v is constant. 

4. THE FRACTION OF DARK HALO LENSES 

According to the model described in §2, we predict the frac- 
tion of dark halo lenses. In our model, there are three parame- 
ters which govern effects of baryonic compression; 6h, M^, and 
7v. Since the constraint from the velocity function of galaxies 
(§3.1) indicates that 5^ should be restricted around ~ 0.78, 
we fix Sh = 0.78 in the remainder of this paper. On the other 
hand, the other parameters, Mj, and were poorly constrained 
from the velocity function of galaxies and the distribution of 
image separations. Therefore we should see the dependence of 
the fraction of dark halo lenses on both Mh and j^. Another 
important parameter we should examine is the inner slope of 
density profile (a in eq. [2]), because it is also known that the 
number of dark halo lenses is extremely sensitive to a (Wyithe 
et al. 2001; Keeton & Madau 2001; Takahashi & Chiba 2001; 
Li & Ostriker 2002; Oguri et al. 2002). 

The fraction of dark halo lenses at image separations 6\ < 
^ < 02 is given by 

f \ - -PdarkO ^l)-Pdark(> Ol) 

where Pdark(> 0) and P(> 0) are calculated from equations (10) 
and (15), respectively. Again, /dark is almost independent of 



the source population (source redshift and flux distribution) 
and cosmological parameters because these primarily change 
the normaUzation of the overall lensing rate which we never 
use. For the range of image separations, we consider following 
two cases: (1) 0.3" < 6 < 15". In this range of image sepa- 
rations, CLASS survey is complete (Helbig 2000; Phillips et 
al. 2001). Among 18 lenses observed in CLASS survey, virtu- 
ally all lenses are known to be galaxy lenses (Rusin 2002). (2) 
0.3" < < 3". Most of current ambiguous quasar pairs which 
may be dark halo lenses have separations 6 > 3" and almost all 
lenses with smaller separations are known to be produced by 
normal galaxies (Kochaneketal. 1999; Rusin 2002). Therefore, 
in this case we can use all lenses with separations 0.3" <6 <3" 
in comparison with our result. Below we assume that there is 
no dark halo lens for both cases. 

Figure 6 plots contours of /dark in Mt^-Jv plane. The con- 
straint from the velocity function of galaxies is also shown. 
From this figure, it is found that the fraction of dark halo lenses 
is sensitive to 7^: /dark increases as 7,, decreases. The reason 
is that the cross section of SIS lenses scales as ctsis oc 7^ (see 
eq. [11]) and the number of galaxy lenses decreases as 7v de- 
creases. The fraction of dark halo lenses also increases as Mh 
decreases, because the number of dark halos which act as dark 
halo lenses becomes large. Figure 6 indicates that /dark cannot 
be determined uniquely, even if we restrict Mh and 7v such that 
they satisfy the constraint from the velocity function of galax- 
ies. This fact is clearly shown in Figure 7. In Figure 7, we plot 
/dark as a function of ■jy. The parameter Mh is chosen so that 
it satisfies the constraint from the velocity function of galaxies. 
This figure indicates that /dark is sensitive to jy. Even if we re- 
strict 1 < 7v < 2 as discussed in §3.3, /dark has uncertainty of 
about one order of magnitude. Figure 7 also shows that /jark 
is sensitive to a. From a = 1.0 to 1.5, /dark also changes about 
one order of magnitude. Therefore, our conclusion is that ro- 
bust predictions of /dark need information on both jy and a. 
If the density profile of dark halos is determined from several 
observations, we can break the degeneracy between Mh and 7v 
from the observation of /dark- 

Although /dark strongly depends on jy, we obtain the lower 
Umit /dark.min by restricting jy. Discussions in §3.3 suggest that 
it is safe to adopt 7,, < 2. Thus by setting 7^ = 2 we can derive 
/dark.min, which is shown in Figure 8. In this figure /dark,min is 
plotted against the inner slope a. The lack of dark halo lenses 
places upper Hmit on the value of /dark which is also shown in 
Figure 8. From this figure, we find that too cuspy inner slope 
of dark halos (a > 1.5) is inconsistent with the lack of dark 
halo lenses: even if we adopt sufficiently large 7v, 7v = 2, we 
predict too much dark halo lenses to reconcile with the observa- 
tion. This constraint has somewhat different meaning from the 
one derived from the statistics of wide separation lensing (Kee- 
ton & Madau 2001 ; Li & Ostriker 2002) because such statistics 
may suffer from the uncertainties of cosmological parameters 
and the source population. 

Are there any other effects which may change /dark? We con- 
sider following two effects which may change our results. One 
is the effect of "substructures" (i.e., galaxies in groups and clus- 
ters) as discussed in §2.2. Substructures may affect the fraction 
of dark halo lenses directly, because it changes only the num- 
ber of galaxy lenses. Moreover the situation that lens galaxies 
which lie in groups or clusters is not rare (e.g., Q0957H-561), 
thus the effect should be addressed quantitatively. Although 
the effect of substructures seems difficult to estimate, Keeon, 
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Fig. 6. — Contour plots for the fraction of dark halo lenses m Mh-7v plane. Thick lines are 68% contours derived from the velocity function (see Fig. 1). The 
inner cusps of darkhalos (eq. [2]) are a = 1.0 (left) and a = 1.5 (right). For the range of image separations, 0.3" <6 < 15" (upper) and 0.3" < 5 < 3" (lower) are 
considered. 



Christlein, & Zabludoff (2000) calculated the fraction of lens- 
ing galaxies which lie in groups and clusters. They predicted 
that ~ 25% of lens galaxies are hkely to be in groups or clusters. 
Following this, the effect of substructures are simply estimated 
by replace the probability of galaxy lensing as Pgai (4/3)Pgai. 
This decreases /dark by a factor 3/4 at most. The fraction of 
dark halo lenses including this effect is also displayed in Figure 
8. As seen in the figure, our result that a> 1.5 is inconsistent 
with the observation is not so affected by the effect of substruc- 
tures. The other effect is the existence of "empty halos", that 
is, dark halos which do not host a galaxy. We have adopted the 
galaxy formation probability Pg(M) of the form (1), and this 
means that we have assumed there is no empty halo at M < Mh. 
Although this assumption does not conflict with the observa- 
tion of galaxy-galaxy lensing (Seljak 2002b), it is still possible 
that the significant amount of empty halos exists at galaxy-mass 
scale. This effect, however, only increases /dark- Therefore 
lower limit of /dark used in Figure 8 is never changed by this 
effect. 

The main result of Wyithe et al. (2001) and Keeton & Madau 
(2001) is that the lensing probabihty of dark halos depends 
strongly on the choice of concentration parameters as well as 
inner slopes a. Therefore, in Figure 9 we examine /dark.min 
for different concentration parameters. In calculating this, we 
adopt the concentration parameter of the following form: 

ritf \ '-norm / ™ \ /•ir>\ 

c^iM,z)~^[-,,^) , (19) 
instead of equation (3), and regard Cnorm as a free parameter. 



We plot /dark.min in ct-Cnoim plauc. From this figure, it is clearly 
seen that there is a strong degeneracy between a and Cnorm as 
reported by Wyithe et al. (2001) and Keeton & Madau (2001). 
Thus actually we can constrain the combination of a and Cnorm. 
or the core mass fraction proposed by Keeton & Madau (2001), 
instead of a and Cnorm separately. 

5. CONCLUSION 

We have studied the effect of baryonic compression assum- 
ing the simple parametric model used by Kochanek (2001). Our 
model has following two elements: the galaxy formation prob- 
ability Pg(M) (eq. [1]) which describes the global efficiency 
of baryonic compression, and the ratio of circular velocities 
of galaxies to virial velocities of dark halos 7v = Vc/vyir which 
means how the inner structure of dark halos is modified due to 
baryonic compression. The model parameters are constrained 
from the observed velocity function of galaxies and the distri- 
bution of image separations in strong gravitational lensing, al- 
though the strong degeneracy between model parameters still 
remains. By using this model, we predict the fraction of dark 
halo lenses /dark (eq. [18]). Here dark halo lenses mean the lens 
systems which are produce by the gravitational potential of dark 
halos. The fraction of dark halo lenses is independent of the 
normalization of total lensing rate, thus is insensitive to cosmo- 
logical parameters and information of sources such as redshift 
and flux distribution. Instead the fraction of dark halo lenses 
is expected to have information on both the effect of baryonic 
compression and the density profile of dark halos such as the 
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Fig. 7. — The fraction of dark halo lenses as a function of 7v. Each solid line is the result with particular choice of a. For each 7v, the value of Mh is chosen so 
that it satisfies the constraint from the velocity function. Observational upper limits of /dark at la level (all lenses and the CLASS survey lenses) are also shown by 
dotted Unes with downward arrows. 



inner density profile for which we modeled p cx r"" . We found 
that /dark is indeed sensitive to the inner slope a, concentration 
parameter Cnom. and model parameters such as jy. Therefore, 
definite predictions of /dark need correct knowledge of baryonic 
compression as well as the density profile of dark halos. This 
also means that we can constrain the model of baryonic com- 
pression from the observation of /dark if the density profile of 
dark halos is well known. 

Although the fraction of dark halo lenses is difficult to pre- 
dict, we can still derive the lower limit of /dark by restricting 
7v < 2 which is inferred from various theories and observations 
(see §3.3). We found that the steep inner profiles (a > 1.5) or 
too centrally concentrated dark halos are inconsistent with the 
lack of dark halo lenses in observations. As described above, 
our result is quite insensitive to cosmological parameters and 
source population. Therefore, our result is complementary to 
the result of Keeton & Madau (2001) who obtained similar con- 
straint on the density profile of dark halos using the total lensing 
probability which also depends strongly on cosmological pa- 
rameters. One of possible systematic effects which may change 
/dark is the effect of galaxies in groups or clusters. By using the 
result of Keeton et al. (2000), we found that this effect changes 
/dark by a factor 3/4 at most. Therefore this effect does not 
change our main result so much. Other important systematic 
effect is the existence of empty halos which are neglected in 



our model. This effect is, however, not important for the lower 
limit of /dark, because this effect only increases the fraction of 
dark halo lenses. 

One possible criticism of our result is that the baryon com- 
pression model we use in this paper is too simple. Many of 
the previous work, however, have not addressed this problem. 
For example, most work of gravitational lensing statistics which 
use the mass function of dark halos assumed that circular veloc- 
ities of galaxies are the same as virial velocities of dark halos 
(this corresponds to 7i, = 1 in our model). But this assumption 
seems to be invalid (e.g., Seljak 2002a), and since the num- 
ber of galaxy lenses scales as cx 7^ the deviation from 7y = 1 
should not be dismissed. This fact is seen even in our simple 
model where the connection between galaxy lenses and dark 
halo lenses sensitively depends on 7,., and is indeed difficult to 
be determined. The importance of baryonic compression, how- 
ever, means that we can constrain the model of baryonic com- 
pression from observations of the fraction of dark halo lenses. 
Although the current sample of gravitational lensing may be too 
small for this purpose, the larger lens sample obtained by e.g., 
SDSS can strongly constrain the model of baryonic compres- 
sion as well as the density profile of dark halos. 

The author would hke to thank Yasushi Suto, Atsushi Taruya, 
and Mamoru Shimizu for useful discussions and comments. 
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Fig. 8. — The lower limit of /dart as a function of inner slope a, when 7v is restricted to 7v < 2. Dashed lines are calculated simply by setting 7v = 2 in Figure 
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Fig. 9. — Contour plots for the lower limit of /dark in the a-Cnonn plane. The effect of galaxies in groups/clusters is taken into account. Dotted lines show the la 
constraint from the lack of dark halo lenses in observations. Dashed lines indicate the value of Cnorm we used in the previous plots. 
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